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The temperature dependence of magnetization and specific heat of the ferromagnetic compound 
DyFe2Zn2o has been measured in detail in various magnetic fields. We have observed anomalous 
magnetic behavior, i.e., a strong anisotropy at 2 K, disappearance of this anisotropy between ap- 
proximately 30 K and T c , and anomalous behavior of the specific heat in the magnetic fields near 
20 K. These anomalous phenomena have been analyzed based on the strong exchange interaction 
between the Fe itinerant electrons and the Dy localized electrons as well as the crystalline electric 
field, Zeeman energy, and an usual exchange interaction between two Dy atoms. The higher T c of 
DyFe2Zn2o compared with that of DyRu2Zn2o is caused by this exchange interaction between the 
Fe and Dy atoms. 
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I. INTRODUCTION 

The series of cubic RX2Zn2o has been recently closely 
examined, where R is a rare earth atom and X is Fe, 
Co, Ru, etc.flHl]. R is located in the cubic symmetric 
site surrounded by 16 Zn atoms. The exchange interac- 
tion between two R atoms is weak because R atoms are 
diluted in this compound. In fact, in the case of X = 
Co or Ru, the Curie temperature is less than 10 K for 
any R. In the case of the Fe series, however, T c is signif- 
icantly enhanced. It has been suspected that the mag- 
netism of Fe atoms is involved in this high T c . However, 
YFe2Zn2o is a nonmagnetic material, i.e., the Fe atoms 
in YFe2Zn2o are non-magnetic. At present, what causes 
the high T c of RFe2Zn2o is not clear. We closely exam- 
ined the magnetization, the magnetic susceptibility and 
the specific heat of DyFc2Zn 2 o in various magnetic fields, 
and found experimentally unusual magnetic anisotropics, 
i.e., the large magneto-crystalline anisotropy at 2 K and 
a metamagnetic transition in the fields along the [100] 
direction, which disappears in the temperatures above 
30 K below T c . These anomalous features are analyzed 
based on the following hamiltonians; the crystalline elec- 
tric field, the exchange interaction between two R atoms 
and the exchange interaction between Fe and R atoms. 
We have taken into account the contribution of the 3d 
electrons of Fe atoms based on the simplified SCR theory. 
In this paper, we show that the enhancement of T c , the 
strong magnetic anisotropy, the metamagnetic behavior 
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at low temperatures, the disappearance of the magnetic 
anisotropy below T c down to 30 K, and no anomaly of 
magnetic specific heat at T c , etc. can be explained by 
our calculation. 



II. EXPERIMENTAL RESULTS 

Figure Q] shows the magnetization curves M(H) in the 
magnetic fields H along the three principal axes [100], 
[110] and [111] at 2 and 30 K. M(H) at 2 K is signif- 
icantly anisotropic, as has been shown by Jia et a2.Q]]. 
The easy axis of magnetization is the [111] direction. At 
zero field, the values of residual magnetic moments M 
along the [110] and [100] directions approximately equal 
to the projection of M to these directions; the calculated 
ratio of M along the [111], [110] and [100] directions is 
l:-\/2/3:l/v / 3. In the figure, a metamagnetic behavior of 
M is seen in H along the [100] direction in approximately 
1.2 T. The anisotropic behavior of M(H) disappears at 
30 K even though this temperature T is much lower than 
T c . In the inset of Fig. [1] the details of the field depen- 
dence of M{H) are shown in H along the [100] direction. 
The field of the metamagnetic transition gradually de- 
creases with increasing temperature, and the metamag- 
netic behavior disappears at 25 K. At 15 and 20 K, the 
metamagnetic behavior is not clearly observed, but an 
unusual M{H) remains in H lower than 0.2 T. 

Figure [5] shows the temperature dependence of 
M(T)/H in H II [100], [110], and [111]. H is applied 
at 0.1 T. From this figure, it is apparent that the strong 
anisotropy of M{T) at low T rapidly shrinks at approxi- 
mately 20 K and disappears above 30 K. The Curie tem- 
perature T c is 53 K, which was determined by M(T) 
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FIG. 1: Magnetization curves M(H) of DyFe 2 Zn 2 o at 2 K 
and 30 K in H || [100], [110], and [111]. Inset shows the 
temperature dependence of A'I(H) in H | [100]. 
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FIG. 2: Temperature dependence of M/H of DyFe2Zn2o in 
H || [100], [110], and [111]. Inset shows the temperature de- 
pendence of H/M. 



curve in H — 0.01 T, because the temperature change 
of M(T) at T c is more clearly observed in low H than 
that in higher H. The inset of Fig. 2 shows the temper- 
ature dependence of H/M, which obeys the Curie- Weiss 
law. The paramagnetic Curie temperature p is 21.1 K, 
and the effective moment /z e ff is 10.8 /ie, which is close 
to the effective moment of the free ion of Dy atom, 10.6 
/ib- The magnetism of DyFe2Zn2o seems to be originated 
from the magnetic ions of rare earth Dy atoms. However 
the temperature dependence of M(T) is unusual; M de- 
creases linearly with increasing T, and the temperature 
change of M at T c is not clear. 

The specific heat C of DyFe2Zn2o and YFe2Zn2o was 



measured. The magnetic part C m of DyFe2Zn2o was de- 
fined as C m = CDyFc 2 Zn 20 ~ CYFc 2 Zn 20I where any correc- 
tions due to the difference of molar weight have not been 
made. Our C m (T) of DyFe2Zri2o was almost the same as 
the one reported by Jia et a2.[lj. The two characteristic 
features of C m (T), which is not shown in the present pa- 
per, are (1) a broad peak around 15 K, (2) an unusual 
and unclear change of C m at T c . 

The specific heat C m was measured in H along the 
three principal axes, and we have found anomalous be- 
havior of C m especially below 20 K. Figure 3 shows the 
temperature dependence of C m (T) in H along the [100] 
direction. In H = 0, C m (T) shows only a broad peak 
with the amplitude of 13 J/mol K at approximately 12 
K. In 0.2 T, however, a tiny thorn appears at 16 K, and 
then the thorn grows into an apparent peak with increas- 
ing H, as shown in the figure. Then the peak abruptly 
disappears in H = 1.4 T. The entropy S m {T) reaches to 
14.2 J/mol K (R In 5.6) at 20 K in H = 0, and decreases 
to 9.0 J/mol K {R In 3) at 20 K in H = 8 T, where R is 
the gas constant. In the case of H along the [111] direc- 
tion, contrary to C m (T) in H // [100], C m (T) steadily 
decreases without any anomalies by applying H, and the 
broad peak is gradually weakened with increasing H and 
disappears in 8 T. 



III. ANALYSES AND DISCUSSION 

We analyze the anomalous behavior of M(T), M(H) 
and C m (T) of DyFe2Zn2o taking into account the mag- 
netism of both the Dy and Fe atoms. First, we esti- 
mate the contribution from the Fe atoms based on the 
experimental results of YFc2Zn.2o, which is a so-called 
nearly-ferromagnet. Figure 4 and the insets show the 
inverse magnetic susceptibility l/x(T), the magnetiza- 
tion curves M(H) and the specific heat divided by T, 
C(T)/T. Note that M is the magnetization per one mole 
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FIG. 4: Temperature dependence of the inverse magnetic sus- 
ceptibility of YFe2Zn2o in H || [111]- Open circles and the 
solid line denote the experimental data and theoretical line, 
respectively. Lower inset shows the temperature dependence 
of M(H) curves at various T in H || [111]. Solid lines de- 
note the theoretical lines. Upper inset shows the specific heat 
divided by T. The line in this inset is a guided line for the 
eyes. 



of YFe2Zn2o, that is, per two Fe atoms. M(H) in the 
inset in Fig. 4 is experimental raw data, but the \ is 
corrected by subtracting small ferromagnetic impurities 
less than 0.0005 ^e/2Fe. The electronic specific-heat co- 
efficient 7o of YFe2Zn 2 o is 0.05 J/mol K, as shown in the 
inset. We analyzed the magnetic properties of YFe2Zn 2 o 
based on the equation, 
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where Mp c denotes the magnetization of Fe atoms. This 
equation is equivalent to the Stoner model when (to 2 ) = 
0. The parameter / is a Stoner's enhancement factor, 
and (to 2 ) is a correlation function of magnetization of 
Fe. Thus, this equation is equivalent to the SCR the- 
ory if this term is correctly taken into account [9l-[lTj|. I n 
this paper, however, to simplify the analysis, we assume 
that (to ) is approximately proportional to the temper- 
ature T, that is, b{m 2 ) = cT. This assumption can hold 
in the paramagnetic region, because the SCR theory pro- 
vides the same form as the Curie- Weiss law in this region. 
Therefore, experimentally, we can use the equation, 



H = (a + cT + dT 2 )M Fc 



(2) 



where a = l/xo — I- In the limit of H — > 0, M Fc — » 0, and 
thus eq. 2 corresponds to the Curie- Weiss law (CW) of 
the localized model if d = 0. The term of d is a correction 
term from the straight line of the CW law. These param- 
eters are experimentally estimated as; a = 96 mol/emu, 



c = 2.18 mol/emu K, d = 3.7 x 10~ 4 mol/emu K 2 . The 
solid line in Fig. 4 is a calculated line using these pa- 
rameters. As shown in the inset in Fig. 4, at low tem- 
peratures in finite magnetic fields, M(H) curve is not a 
straight line. The second term in eq. 2 is a correction 
term to correct the straight line. Thus, b is experimen- 
tally obtained as 9.19xl0~ 5 Oe/(emu/mol) 3 . The solid 
lines of M(H) curves in the inset of Fig. 4 are the cal- 
culated ones using these parameters. The parameter d is 
not effective at low temperatures. The parameter \o is 
estimated as 6.85 x 10 -4 emu/mol using the experimental 
value of 7o. The Stoner's enhancement factor a = xol is 
0.93. 

Next, we define the following hamiltonian to analyze 
the anomalous properties of DyFe2Zn2o, 



H — HqeF + -ffcxch + Hz + H 



FcR: 



(3) 



where Hqef is the crystalline-electric-field (CEF) hamil- 
tonain, Hz the Zeeman hamiltonian, and i? cxc h the ex- 
change hamiltonian between the two Dy atoms. The 
former three terms are hamiltonians for Dy atoms in 

DyFe 2 Zn 2 o. ^cef is written as VF(^0 4 + T^Oa), H z 
is -M n H ext , and H cxch = -M R H moh where M R is 
the magnetization of Dy atoms, H oxt the external field, 
H mo \ the molecular field at Dy atom caused by surround- 
ing Dy atoms, M R = -gjfi B J, H mol = n RR M R , and 
n RR the exchange constant between two Dy atoms. The 
last term in eq. 3 is an exchange interaction between Dy 
and Fe atoms, which is written as 
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where M Fc is the magnetization of Fe atoms, and n Fe R 
the exchange constant between Dy and Fe atoms. Note 
that the physical quantities of Af R , J, H mo i, and H CKt 
are three-dimensional vectors. The molecular field at the 
Dy atom caused by the surrounding Fe atoms is written 
as 



^FcrMfo = n FoR Af Fc 



H, 



n FeR M R 



71F c rM R | 



(5) 



where M Fc — |Mf c |, and the term in the parenthesis 
is a unit vector parallel to M Fc . In eq. 5, we assumed 
that the magnetic moment of Fe is isotropic, i.e., aligning 
parallel to the direction of H ext +71f rM r , which is the 
sum of Hext and the molecular field at the Fe site. The 
direction of H cxt is not always parallel to the direction of 
M R , especially in the ferromagnetic temperature region. 
Consequently, M Fc obeys the following equation, 



tifcrMrI = {a + cT + dT 2 )M Fc + 6M| C .(6) 



TABLE I: Characteristic parameters. a [mol/emu], c 
[mol/emu K], d [lCT 4 mol/emu K 2 ], b [10" 5 Oe/(emu/mol) 3 ], 
9 P [K], Peff [Ms/Fe], a: Stoner's enhancement factor. 
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b a 


Op Pcff 


#2 96 2.18 


-3.7 


9.19 0.93 


-44.1 1.36 
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We summarize the hamiltonian for Dy atom in the 
frame of the molecular-field approximation as follows: 



H = i?cEF + 9jMb(1 + h)H cxt J 

-(ffJA*B) 2 (nRR + fexch){J)J, 



where 



h = 

fcxch — 



™FcrMf c 



|-ffext - gjflBnFeR(J) 
n FeB. M Fc 



The magnetization of Dy atom is calculated by 

Tr Jcxp(-/3g) 
[ ' Tr exp(-/3/f) ' 

where /3 — l/k^T. The specific heat C is calculated by 



(7) 

(8) 
(9) 

(10) 



(ii) 



where the second term in the parenthesis is a correction 
term to correct the double-counting of the exchange en- 
ergy between the two Dy atoms. 

There are eight parameters to fit the experimental 
M(T), M(H), C(T), and C(H) in this calculation: 
x, W, nRR, np C R,a, c, d and b. However, the param- 
eters a, c, d and b are fixed as in YFe2Zii2o- More- 
over, the parameters x, W and tirr are estimated by 
adjusting the magnetic crystalline anisotropy and T c for 
DyRu2Zn2o, which is the compound without magnetic 
Fe atoms, where T c of DyRu2Zn2o is 4 K. The final pa- 
rameter is npoR, which is determined by adjusting T c of 
DyFe2Zn2o- The eight fitting parameters seem to be eas- 
ily found. 

However, difficulties of the calculation exist in the it- 
eration procedures when we solve simultaneously both 
Mr of Dy atoms and Mp c of Fe atoms using eqs. 6 and 
10. In this paper, we have adopted a two-step procedure 
to solve these equations. First, we ignore the anisotropy 
of Mr. When Hcef is excluded from the total hamilto- 
nian, the directions of Mp c and .Mr are parallel to if ex t 
in the case of positive np R. Mr is easily expressed using 
the Brillouin function Bj(x) as 



(J z ) = -JBj(x), 



where 



gjH B J(H + n RR (J z ) + n FcR M Fc ) 



(12) 



(13) 



In the case of negative tifoR, the directions of Mp and 
Mr are not always parallel to ff ex t- However, we used 
cq. 12 for both cases for simplicity. Using eqs. 6 and 
12, Mpo and (J z ) are easily calculated. As the results 
of this step, the contribution of Mp c is found to be ap- 
proximately 4 % of Mr. Then, as the second step of the 
calculation, the anisotropic J of Dy atom is re-calculated 
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FIG. 5: Calculated M(H) curves at 2 and 30 K in H \\ [100], 
[110], and [111]. Inset shows the calculated M(H) curves at 
various T in H \\ [100]. 
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FIG. 6: Calculated temperature dependence of M/H in H \\ 
[100], [110] and [111]. Inset shows the calculated H/M in H \\ 
[100]. The thin line in the inset is a guided line for the eyes. 



based on eq. 10 by the iteration procedure. In this step, 
we fixed Mp G at the value calculated in the first step. 

Figures 5, 6 and 7 show the calculated results of 
M R (H), M R (T) and C(T), respectively. The calculated 
M in these figures is Mr, because Mp c is presumably 
smaller than M R . The CEF parameters and the exchange 
interaction constants are listed in Table 2. The energy 
splitting by CEF is considerably small; the ground state 
is Tg, — —7.5 K, the first excited state is Tg, = —3.5 K, 
and the overall energy splitting is 16.3 K. 

The calculated results of M(H) and M(T) are in good 
agreement with the experimental results in Figs. 1 and 
2 both qualitatively and quantitatively. The calculated 
results represent the following characteristic features of 
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FIG. 7: Calculated specific heat C m (T) in various H || [100]. 



increasing H and then disappears above 3 T. In the cases 
of H along the [110] and [111] directions, we could also 
analyze the experimental results of M and C by our cal- 
culation. 

We summarize the experimental and calculated results 
of DyFe2Zn2o- DyFe2Zii2o is a ferrimagnet with two mag- 
netic moments of Dy and Fe atoms. The former ion has 
10 /ie/Dy, the latter 0.4 /i B /2Fe. This small magnetic 
moment of Fe atom has been currently inferred from the 
Mossbauer experiments [l2j . The overall splitting of CEF 
is 16 K. The exchange interaction tirr between the two 
Dy ions is extremely weak compared with the exchange 
interaction ?if c r between the Dy and Fe atoms. Thus 
the T c of DyFe2Zn2o is enhanced by this riFeR, because it 
is 80 times larger in magnitude than tirr in DyFc2Zii2o- 
This small urr is caused by the distance between Dy ions 
in DyFe2Zn2o, and the large value of n^eR is presumably 



DyFe2Zn2o: the calculated anisotropy M(H) at 2 K coin- 
cides with the experimental one, the easy axis of magne- 
tization is parallel to [111], the metamagnetic transition 
occurs in H along the [100] direction in approximately 
2 T at 2 K, the metamagnetic transition field decreases 
with increasing T and goes down to zero near 20 K, the 
magnetic anisotropy disappears for all the principal axes 
above 30 K, Concerning C(T) in H along the [100] di- 
rection, the calculated results reproduce well the exper- 
imental features: C(T) in H = makes a broad peak 
at approximately 10 K without any anomalous thorn, a 
small anomalous thorn appears by applying H along the 
[100] direction and it grows into an apparent peak with 



TABLE II: CEF parameters, W [K] and x, and the ex- 
change constants between R-R and R-Fe, urr [T//j,b} and 
"■Fell [T/^b]- 

W x TtRR n Fc R 

-0.036 -0.2 0.188 -14.8 

due to the large overlapping of the wave function between 
the itinerant 3d and the localized 4f electrons. 
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